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Abstract In the classical (min-cost) Steiner tree problem, we are given 
an edge- weighted undirected graph and a set of terminal nodes. The goal 
is to compute a min-cost tree 5* which spans all terminals. In this pa- 
per we consider the min-power version of the problem (a.k.a. symmetric 
multicast), which is better suited for wireless applications. Here, the goal 
is to minimize the total power consumption of nodes, where the power 
of a node v is the maximum cost of any edge of 5* incident to v. Intu- 
itively, nodes are antennas (part of which are terminals that we need to 
connect) and edge costs define the power to connect their endpoints via 
bidirectional links (so as to support protocols with ack messages). Ob- 
serve that we do not require that edge costs reflect Euclidean distances 
between nodes: this way we can model obstacles, limited transmitting 
power, non-omnidirectional antennas etc. Differently from its min-cost 
counterpart, min-power Steiner tree is NP-hard even in the spanning tree 
case (a.k.a. symmetric connectivity), i.e. when all nodes are terminals. 
Since the power of any tree is within once and twice its cost, comput- 
ing a Pat < ln(4) -I- e [Byrka et al.'lO] approximate min-cost Steiner tree 
provides a 2pst < 2.78 approximation for the problem. For min-power 
spanning tree the same approach provides a 2 approximation, which was 
improved to 5/3 + e with a non-trivial approach in [Althaus et al.'06]. 

In this paper we present an improved approximation algorithm for 
min-power Steiner tree. Our result is based on two main ingredients. We 
present the first decomposition theorem for min-power Steiner tree, in 
the spirit of analogous structural results for min-cost Steiner tree and 
min-power spanning tree. Based on this theorem, we define a proper 
LP relaxation, that we exploit within the iterative randomized round- 
ing framework in [Byrka et al.'lO]. A careful analysis of the decrease of 
the power of nodes at each iteration provides a31n4— |-|-£ < 1.91 
approximation factor. The same approach gives an improved 1.5 -I- e 
approximation for min-power spanning tree as well. This matches the 
approximation factor in [Nutov and Yaroshevitch'09] for the special case 
of min-power spanning tree with edge weights in {0, 1}. 

1 Introduction 



Consider the following basic problem in wireless network design. We are given a 
set of antennas, and we have to assign the transmitting power of each antenna. 
Two antennas can exchange messages (directly) if they are within the transmis- 
sion range of each other (this models protocols with ack messages) . The goal is 
to find a minimum total power assignment so that a given subset of antennas 
can communicate with each other (using a multi-hop protocol). 
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We can formulate the above scenario as a min-power Steiner tree problem 
(a.k.a. symmetric multicast). Here we are given an undirected graph G — (V, E), 
with edge costs c : E ^ Q>Oi and a subset R of terminal nodes. The goal is to 
compute a Steiner tree S spanning i?, of minimum power p{S) := Xlt;GV(s) Ps{v), 
with psiv) ■= niax„„g^(5){c(uu)}. In words, the power of a node v with respect 
to tree S is the largest cost of any edge of S incident to v, and the power of S 
is the sum of the powers of its nodeqj. The min-power spanning tree problem 
(a.k.a. symmetric connectivity) is the special case of min-power Steiner tree 
where R = V, i.e. all nodes are terminals. Let us remark that, differently from 
part of the literature on related topics, we do not require that edge costs reflect 
Euclidean distances between nodes. This way, we are able to model obstacles, 
limited transmitting power, antennas which are not omnidirectional, etc. 

The following simple approximation-preserving reduction shows that min- 
power Steiner tree is at least as hard to approximate as its min-cost counterpart: 
given a min-cost Steiner tree instance, replace each edge e with a path of 3 edges, 
where the boundary edges have cost zero and the middle one has cost c(e)/2. 
Hence the best we can hope for in polynomial time is a c approximation for 
some constant c > 1. It is known [llllj that, for any tree S of cost c{S) := 
See-E(5') c(e), 

ciS)<p{S)<2c{S). (1) 

As a consequence, a pst approximation for min-cost Steiner tree implies a 2pst 
approximation for min-power Steiner tree. In particular, the recent ln(4) -|- £ < 
1.39 approximatiorQ in [3| for the first problem, implies a 2.78 approximation 
for the second one: no better approximation algorithm is known to the best of 
our knowledge. 

Differently from its min-cost version, min-power spanning tree is NP-hard 
(even in quite restricted subcases) [1111) . By the above argument, a min-cost 
spanning tree is a 2 approximation. However, in this case non-trivial algorithms 
are known. A l-|-ln2-|-£ < 1.69 approximation is given in [5] . This was improved 
to I -I- £ in [I]. If edge costs are either or 1, the approximation factor can be 
further improved to | -I- £ [12] . Indeed, the same factor can be achieved if edge 
costs are either a or b, with < a < 6: this models nodes with two power states, 
low and high. All these results exploit the notion of A:-decomposition. The first 
result is obtained with a greedy algorithm, while the latter two use (as a black 
box) the FPTAS in |13j for the min-cost connected spanning hypergraph problem 
in 3-hypergraphs. We will also use /c-decompositions, but our algorithms are 
rather different (in particular, they are LP-based). 

Our Results. In this paper we present an improved approximation algorithm 
for min-power Steiner tree. 

Theorem 1. There is an expected 3\n 4 — j+e < 1.909 approximation algorithm 
for min-power Steiner tree. 

^ When S is clear from the context, we will simply write p{v). 
^ Throughout this paper e denotes a small positive constant. 
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Our result is based on two main ingredients. Informally, a /c-decomposition of 
a Steiner tree is a collection of (possibly overlapping) subtrees of S, each 
one containing at most k terminals, which together span S. The power/cost 
of a decomposition is the sum of the powers/costs of its component^. It is 
a well-known fact (see [2 and references therein) that, for a constant k large 
enough, there exists a fc-decomposition of cost at most (1 + e) times the cost of 
S. A similar result holds for min-power spanning tree [T]. The first ingredient in 
our approximation algorithm is a similar decomposition theorem for min-power 
Steiner tree, which might be of independent interest. This extends the qualitative 
results in |ll2j since min-power Steiner tree generalizes the other two problems. 
However, the dependence between e and k is worse in our construction. 

Theorem 2. (Decomposition) For any h > 3 and any Steiner tree S, there 
exists a -decomposition of S of power at most (1 + ^)p{S). 

Based on this theorem, we are able to compute a 1 -I- e approximate solution 
for a proper component-based LP-relaxation for the problem. We exploit this 
relaxation within the iterative randomized rounding algorithmic framework in 
[3]: we sample one component with probability proportional to its fractional 
value, set the corresponding edge costs to zero and iterate until there exists a 
Steiner tree of power zero. The solution is given by the sampled components 
plus a subset of edges of cost zero in the original graph. A careful analysis of 
the decrease of node powers at each iteration provides a 31n4 — | -I- e < 1.91 
approximation. We remark that, to the best of our knowledge, this is the only 
other known application of iterative randomized rounding to a natural problem. 

The same basic approach also provides an improved approximation for min- 
power spanning tree. 

Theorem 3. There is an expected |+£ approximation algorithm for min-power 
spanning tree. 

This improves on [13], and matches the approximation factor achieved in [12] 
(with a drastically different approach!) for the special case of 0-1 edge costs. 

Preliminaries and Related Work. Min-power problems are well studied in 
the literature on wireless applications. Very often here one makes the assumption 
that nodes are points in or M"^, and that edge costs reflect the Euclidean 
distance d between pairs of nodes, possibly according to some power law (i.e., 
the cost of the edge is d'^ for some constant c typically between 2 and 4). This 
assumption is often not realistic for several reasons. First of all, due to obstacles, 
connecting in a direct way geographically closer nodes might be more expensive 
(or impossible). Second, the power of a given antenna might be upper bounded 
(or even lower bounded) for technological reasons. Third, antennas might not 
be omnidirectional. All these scenarios are captured by the undirected graph 
model that we consider in this paper. A relevant special case of the undirected 

^ Due to edge duplication, the cost of the decomposition can be larger than c{S). Its 
power can be larger than p{S) even for edge disjoints components. 
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graph model is obtained by assuming that there are only two edge costs a and 
6, < a < 6. This captures the practically relevant case that each node has only 
two power states, low and high. A typical goal is to satisfy a given connectivity 
requirement at minimum total power, as we assume in this paper. However, it 
makes sense also to consider the min-max version of the problem, where one 
wants to minimize the maximum power. 

Several results are known in the asymmetric case, where a unidirectional link 
is established from u to w iff v is within the transmission range of u (and possibly 
the vice versa does not hold). For example in the asymmetric unicast problem 
one wants to compute a min-power directed path from node s to node t. This 
problem can be solved in polynomial time, say, via dynamic programming. In the 
asymmetric connectivity problem one wants to compute a min-power spanning 
arborescence rooted at a given root node r. This problem is NP-hard even in 
the 2-dimensional Euclidean case [7], and a minimum spanning tree provides a 
12 approximation for the Euclidean case (while the general case is log-hard to 
approximate) [14| . The asymmetric multicast problem is the generalization of 
asymmetric connectivity where one wants a min-power arborescence rooted at r 
which contains a given set R of terminals. As observed in [T], the same approach 
as in [14j provides a 12pat approximation for the Euclidean case, where pst is 
the best-known approximation for Steiner tree in graphs. In the complete range 
assignment problem one wants to establish a strongly connected spanning sub- 
graph. The authors of [11 present a 2-approximation which works for the undi- 
rected graph model, and show that the problem is NP-hard in the 3-dimensional 
Euclidean case. The NP-hardness proof was extended to 2 dimensions in [5]. 
Recently, the approximation factor was improved to 2 — 5 for a small constant 
(5 > with a highly non-trivial approach [1]. The same paper presents a 1.61 
approximation for edge costs in {a, 6}, improving on the 9/5 factor in |6]. 

In this paper we consider the symmetric case, where links must be bidirec- 
tional (i.e. u and v are not adjacent if one of the two is not able to reach the 
other). This is used to model protocols with ack messages. The symmetric uni- 
cast problem can be solved by applying Dijkstra's algorithm to a proper auxiliary 
graph [1] . The symmetric connectivity and multicast problems are equivalent to 
min-power spanning tree and min-power Steiner tree, respectively. We already 
discussed the known results on these problems. One can also consider higher 
connectivity requirements. For example, 0{k) and 0(log^ n) |10j approxima- 
tion algorithms are known for the problem of computing a min-power /c- vertex 
connected subgraph (with bidirectional links). 

Proofs which are omitted due to lack of space are given in the appendix. The 
min-power Steiner tree is denoted by S* . 

2 A Decomposition Theorem for Min-Power Steiner Tree 

A k- component is a tree which contains at most k terminals. If internal nodes 
are non-terminals, the component is full. A k- decomposition of a Steiner tree 
S over terminals i? is a collection of ^-components on the edges of S which 
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span S and such that the following auxiliary component graph is a tree: replace 
each component C with a star, where the leaves are the terminals of C and the 
central node is a distinct, dummy non-terminal vc- Observe that, even if the 
component graph is a tree, the actual components might share edges. When the 
value of k is irrelevant, we simply use the terms component and decomposition. 
We will consider fc-decompositions with k = 0(1). This is useful since a min- 
power component C on a constant number of terminals can be computed in 
polynomial tim^. The assumption on the component graph is more technical, 
and it will be clearer later. Intuitively, when we compute a min-power component 
on a subset of terminals, we do not have full control on the internal structure of 
the component. For this reason, the connectivity requirements must be satisfied 
independently from that structure. 

Assume w.l.o.g. that S consists of one full component. This can be enforced 
by appending to each terminal v a dummy node v' with a dummy edge of cost 
0, and replacing v with v' in the set of terminals. Any decomposition into (full) 
/c-components of the resulting tree can be turned into a fc-decomposition of the 
same power for S by contracting dummy edges, and vice versa. 

Next lemma shows that one can assume that the maximum degree of the 
components in a decomposition can be upper bounded by a constant while losing 
a small factor in the approximation (see also Figure [1} . 

Lemma 1. For any A>3, there exists a decomposition of S of power at most 
(1 + p^/2]-i )p(^) whose components have degree at most A. 

Proof. The rough idea is to split S at some node v not satisfying the degree 
constraint, so that the duplicated copies of v in each obtained component have 
degree (less than) A. Then we add a few paths between components, so that the 
component graph remains a tree. All the components but one will satisfy the 
degree constraint: we iterate the process on the latter component. 

In more detail, choose any leaf node r as a root. The decomposition initially 
consists of S only. We maintain the invariant that all the components but possibly 
the component Cr containing r have degree at most A. Assume that Cr has 
degree larger than A (otherwise, we are done). Consider any split node v of 
degree d{v) — d + 1 > A + 1 such that all its descendants have degree at 
most A. Let ui, . . . ,Ud be the children of v, in increasing order of c{vui). Define 
A' e [2,Z\ — 1]. We partition the u^'s by iteratively removing the 

first A' children, until there are at most A — 2 children left: let Vi, . . . , V/i be the 
resulting subsets of children. In particular, for i < h, Vi = . . . , UiA'}. 

For i = l,...,/i — 1, we let Ci be a new component induced by v, Vi, and the 
descendants of Vi. The new root component Ch is obtained by removing from Cr 

* One can guess (by exhaustive enumeration) the non-terminal nodes of degree at least 
3 in C, and the structure of the tree where non-terminals of degree 2 are contracted. 
Each edge vu of the contracted tree corresponds to a path P whose internal nodes 
are non-terminals of degree 2: after guessing the boundary edges uu' and v'v of P 
(which might affect the power of u and v, respectively), the rest of P is w.l.o.g. a 
min-power path between u' and v' (which can be computed in polynomial time [l]). 
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the nodes Ui</iV^(Ci) — {v} and the correspondmg edges. In order to mamtain 
the connectivity of the component graph (which might be lost at this point), 
we expand C^+i, i > 1, as foUows: let Pj be any path from v to some leaf 
which starts with edge vuj. We append to Ci+i the path Pm{i) which minimizes 
p{Pj) — c{vuj) over j G Vi. After this step, the component graph is a tree. The 
invariant is maintained: in fact, the degree of any node other than v can only 
decrease. In each Ci, i < h, v has degree either Z\' or Z\' + 1, which is within 2 
and A. Since A-2-A' < \Vh\ < A-2, the cardinality |V,J+2 of u in is also 
in [2,Z\]. By the choice of v, all the components but Cr have maximum degree 
A. Observe that Cr loses at least A' — 1 > 1 nodes, hence the process halts. 

In order to bound the power of the final decomposition, we use the following 
charging argument, consisting of two charging rules. When we split Cr at a given 
node u, the power of v remains p{v) in Ch and becomes c{vuiA') in the other 
Ci's. We evenly charge the extra power c{vuiA') to nodes Vi+i: observe that each 
u, e is charged by < < < 

Furthermore, we have an extra increase of the power by p{Pm{i)) ~ c{vum{i)) 
for every i < h: this is charged to the nodes of the paths Pj — {v} with Uj G 

~ {''^m(i)}: in such a way that no node w is charged by more than ^}_i P{w). 
This is possible since there are A' ~ 1 such paths, and the nodes of each such 
path have total power at least p{Pm{i)) — c{vum(i)) by construction. 

Each node w can be charged with the second charging rule at most once, 
since when this happens w is removed from Cr and not considered any longer. 
When w is charged with the first charging rule, it must be a child of some split 
node V. Since no node is a split node more than once, also in this case we charge 
w at most once. Altogether, each node v is charged by at most ^^_^ p{v). □ 

Proof. (Theorem\^ Apply Lemma [1] with A ~ h to hence obtaining a de- 
composition of power at most [|^/2]-i P('^) whose components have degree at 
most h. We describe an decomposition of each such component C with more 
than terminals (see also Figure [5]). Root C at any non-terminal r, and short- 
cut internal nodes (other than r) of degree 2. For any internal node v of C, let 
P{v) be the path from v to its rightmost child r^, and then from to some 
leaf terminal £(w) using the leftmost possible path. Observe that paths P{v) are 
edge disjoint. Pick a value q G {0, 1, . . . , /i — 1} uniformly at random, and mark 
the nodes at level ^ = q (mod h) . Consider the partition of C into edge-disjoint 
subtrees T which is induced by the marked levels. Finally, for each such subtree 
T, we append to each leaf ^; of T the path P{v): this defines a component Cy. 

Trees T have at most leaves: hence components Ct contain at most 
terminals each. Observe that the component graph remains a tree. In order to 
bound the power of components Ct, note that each node u in the original tree 
has in each component a power not larger than the original power p{u): hence it 
is sufficient to bound the expected number p,u of components a node u belongs 
to. Suppose u is contracted or a leaf node. Then u is contained in precisely the 
same components as some edge e. This edge belongs deterministically to one 
subtree T (hence to Ct), and possibly to another component Ct' if the node v 
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with e G Piv) is marked; the latter event happens with probabihty 1/h. Hence 
in this case /i„ < 1 + 1/h. For each other node u, observe that u belongs to 
one subtree T if it is not marked, and to at most two such subtrees otherwise. 
Furthermore, it might belong to one extra component Ct' if the node v with 
ulu € P{v) is marked, where Z„ is the leftmost child of u. Hence, < 1 + 2//i in 
this case. Altogether, the decomposition of C has power at most (1 + 2/h)p{C) 
in expectation. 

From the above discussion, there exists (deterministically) an /i'* decomposi- 
tion of power at most + ji)p{S) < (1 + ^)p{S). □ 

We remark that for both min-cost Steiner tree and min-power spanning tree 
(which are special cases of min-power Steiner tree), improved (l-l-^^^)-approximate 
c'' decompositions, c = 0(1), are known [112] . Finding a similar result for min- 
power Steiner tree, if possible, is an interesting open problem in our opinion (even 
if it would not directly imply any improvement of our approximation factor) . 

3 An Iterative Randomized Rounding Algorithm 

In this section we present an improved approximation algorithm for min-power 
Steiner tree. Our approach is highly indebted to |3,. We consider the following 
LP relaxation for the problem: 

min ^ PQ ■ XQ,s {LPpow) 

(Q,s):a£QCR 

s.t. >1, V0 7^ C i? - {r}; 

(Q,3):aeQCR, 

XQ,s > 0, Vs G Q C 7?. 

Here r is an arbitrary root terminal. There is a variable xq^s for each subset 
of terminals Q and for each s ^ Q: the associated coefficient pq is the power 
of a min-power component Cq on terminals Q. In particular, S* — Cr induces 
a feasible integral solution (where the only non-zero variable is xn^r = !)■ Let 
Cq.s be the directed component which is obtained by directing the edges of Cq 
towards s. For a fractional solution x, let us define a directed capacity reservation 
by considering each {Q,s), and increasing by xq^s the capacity of the edges in 
Cq^s Then the cut constraints ensure that each terminal is able to send one 
unit of (splittable) flow to the root without exceeding the mentioned capacity 
reservation. We remark that the authors of ;_3, consider essentially the same 
LP, the main difference being that pq is replaced by the cost cg of a min-cost 
component on terminals C^. In particular, the set of constraints in their LP is 

® Another technical difference w.r.t. [3] is that they consider only full components: this 
has no substantial impact on their analysis, and allows us to address the Steiner and 
spanning tree cases in a unified way. 
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Algorithm 1 An iterative randomized rounding approximation algorithm for 
min-power Steiner tree. 



(1) For t = 1,2, . . . 

(la) Compute a 1+e approximate solution a;* to LPpom (w.r.t. the current instance), 
(lb) Sample one component C', where C' — Cq with probability 

SsgQ ^Q,s/ X^(Q' s') s^Q' s' . Set to zero the cost of the edges in C' and 

update LPpow 

(Ic) If there exists a Steiner tree of power zero, return it and halt. 



the same as in LPpow This allows us to reuse part of their results and techniques, 
which rely only on the properties of the set of constraint^. 

Given Theorem [21 the proof of the following lemma follows along the same 
line as in [Sj. 

Lemma 2. For any constant e > 0, a 1 + e approximate solution to LPpow can 
be computed in polynomial time. 

We exploit LPpow within the iterative randomized rounding framework in 
[3]. Our algorithm (see also Algorithm [1]) consists of a set of iterations. At each 
iteration t we compute a (l + e)-approximate solution to LPpow, and then sample 
one component C* = Cq with probability proportional to J2si£Q s- We set 
to zero the cost of the edges of C* in the graph, updating LPpow consequentljl^- 
The algorithm halts when there exists a Steiner tree of cost (and power) zero: 
this halting condition can be checked in polynomial time. 

Lemma 3. Algorithmll\ halts in a polynomial number of rounds in expectation. 

4 An Improved Approximation. 

In this section we bound the approximation factor of Algorithm [TJ both in the 
general and in the spanning tree case. Following pj, in order to simplify the 
analysis let us consider the following variant of the algorithm. We introduce a 
dummy variable Xr,r with pr ^ corresponding to a dummy component con- 
taining the root only, and fix X f so that the sum of the x's is some fixed value 
M in all the iterations. For M = n'-^^^) large enough, this has no impact on the 
power of the solution nor on the behaviour of the algorithm (since sampling the 
dummy component has no effect). Furthermore, we let the algorithm run forever 
(at some point it will always sample components of power zero). 

® Incidentally, this observation might be used to address also other variants of the 

Steiner tree problem, with different objective functions. 
^ In the original algorithm in 3 , the authors contract components rather than setting 

to zero the cost of their edges. Our variant has no substantial impact on their analysis, 

but it is crucial for us since contracting one edge (even if it has cost zero) can decrease 

the power of the solution. 
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Let 5* be the min-power Steiner tree at the beginning of iteration t (in par- 
ticular, = S*). For a given sampled component C*, we let p(C*) be its power 
in the considered iteration. We define similarly p{S*) and the corresponding cost 
c{S*). The expected approximation factor of the algorithm is bounded by: 

Hence, it is sufficient to provide a good upper bound on E[p{S^)]. We exploit 
the following high-level (ideal) procedure. We start from S — S* , and at each 
iteration t we add the sampled component C* to S and delete some bridge edges 

in E{S) n E{S*) in order to remove cycles (while maintaining terminal con- 
nectivity). By construction, is a feasible Steiner tree at any time. Furthermore, 
the power of S at the beginning of iteration t is equal to the power p(C/*) of the 
forest of non-deleted edges U* of S* at the beginning of the same iteratiorH. In 
particular, p{S^) < p{U*) = J^vPu'i'")- 

At this point our analysis deviates (and gets slightly more involved) w.r.t. 
[3]: in that paper the authors study the expected number of iterations before 
a given (single) edge is deleted. We rather need to study the behavior of col- 
lections of edges incident to a given node v. In more detail, let el, ... , et^"^ be 

the edges of S* incident to v, in decreasing order of cost c], > >,..., > ct^^^ 
(breaking ties arbitrarily). Observe that Pu* {v) — Cu during the iterations when 
all edges e^, . . . , ej,~^ are deleted and e\, is still non-deleted. Define 5^ as the ex- 
pected number of iterations before all edg deleted. For notational 

convenience, define also 5^ = c^J^^''^^ = 0. Then 

d{v) d(v) 

E[^Put{-v)] = ^cl{Sl -5r^) = ^Sl{cv ~Cv^^)- (3) 

i i = l i = l 

We will provide a feasible upper bound (5* on SI for all v (for a proper choice of 
the bridge edges B*) with the following two properties for all i: 

In words, the (5*'s are increasing (which is intuitive since one considers larger 
sets of edges) but at decreasing speed. Consequently, from ^ one obtains 

d{v) d(v) 

ElJ^Pu^ (y)] < + max{5' - S'-'} V cl = S'cl + {S^ ~ S') V 4. (4) 

— z>2 — — 

i ^ 1=2 i=2 

Inspired by (jlj , we introduce the following classification of the edges of S* . We 
say that the power of node v is defined by ej. We partition the edges of S* into 
the heavy edges H which define the power of both their endpoints, the middle 
edges M which define the power of exactly one endpoint, and the remaining light 
edges L which do not define the power of any node. Let c{H) ~ c{S*) and 
c(M) = 7M c{S*). Observe that ^(5**) = ac{S*) where a = 2-fH + 1m S [1, 2]. 

* Since edge weights of sampled components are set to zero, any bridge edge can be 
replaced by a path of zero cost edges which provides the same connectivity. 



9 



Note also that in Q heavy edges appear twice with coefficient S^, middle edges 
appear once with coefficient 5^ and once with coefficient 6^—6^, and hght edges 
appear twice with coefficient 6"^ — 5^. Therefore one obtains 

t V t 

= {25\h + S^JM + 2(5' - 5')(1 -jH- 7m)) ■ c(S*) 

= (2(5' -S'') + [23' - 5^)a) ■ p{S*). (5) 

Summarizing the above discussion, the approximation factor of the algorithm 
can be bounded by 



MpiS") ^ '"^^ " - Mp{S*) ^ ''^^ " - M 

We next provide the mentioned bounds satisfying Properties (a) and (b): 
we start with the spanning tree case and then move to the more complex and 
technical general case. 



4.1 The Spanning Tree Case. 

Observe that in this case the optimal solution T* :~ S* is by definition a termi- 
nal spanning tree (i.e. a Steiner tree without Steiner nodes). Therefore we can 
directly exploit the following claim in [3]. 

Lemma 4. /5/ Let T* be any terminal Steiner tree. Set T := T* and consider 
the following process. Fort = 1, 2, . . ..• (a) Take any feasible solution x* to LPpowl 
(b) Sample one component C* = Cq with probability proportional to variables 
Xq (c) Delete a subset of bridge edges B* from E{T) n E{T*) so that all the 
terminals remain connected in T — U . There exists a randomized procedure 
to choose the B* 's so that any W C E{T*) is deleted after M H\\y\ iterations in 
expectatioT^. 

By Lemma 2] with W = {e^, . . . , e* }, we can choose 6^ = M ■ Hi. Observe 
that these (5"s satisfy Properties (a) and (b) since ^ ^7*^ ~ ITT ^ positive 
decreasing function of i. Theorem [3] immediately follows by ([6]) since = 

(l+e)MH2 _ 3 
M — \^ ~'~ ' 2- 



4.2 The General Case. 

Here we cannot directly apply Lemma 2] since S* might not be a terminal span- 
ning tree: w.l.o.g. assume that S* consists of one full component. Following [3], 
we define a proper auxiliary terminal spanning tree T*, the witness tree (see 
also Figure [3]). We turn S* into a rooted binary tree S^^^ as follows: Split one 

® Hq := X]i=i T ^he q-th harmonic number. 
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edge, and root the tree at the newly created node r. Split internal nodes of de- 
gree larger than 3 by introducing dummy nodes and dummy edges of cost zero. 
We make the extra assumption^i. that we perform the latter step so that the 
i most expensive edges incident to a given node appear in the highest possible 
(consecutive) levels of 5^j„. Finally, shortcut internal nodes of degree 2. Tree T* 
is constructed as follows. For each internal node v in S^^^-^ with children u and 
z, mark uniformly at random exactly one of the two edges vu and vz. Given 
two terminals r' and r", add r'r" to T* iff the path between r' and r" in Sl^^ 
contains exactly one marked edge. We associate to each edge /' G E{Sl^^) a 
(non-empty) witness set W{f') of edges of T* as follows: e = uv ^ E{T*) be- 
longs to W{f') iff the path between u and v in Sl^^^ contains /'. There is a 
many-to-one correspondence from each / G E{S*) to some /' e E{S^^^): we let 
W{f) := Wif). 

We next apply the same deletion procedure as in Lemma|3]to T* . When all the 
edges in W{f) are deleted, we remove / from S*: this process defines the bridge 
edges i?* that we remove from S at any iteration. As shown in [3j , the non-deleted 
edges J7* of S* at the beginning of iteration t plus the components which are 
sampled in the previous iterations induce (deterministically) a feasible Steiner 
tree. Hence also in this case we can exploit the upper bound p{S*) < p{U^) = 
J2vPu*i'^)- ^'^^ define W^{v) := U'j^iW{ei). In particular, in order to delete 
all the edges el, . . . ,e* we need to delete W'^{v) from T*. The next technical 
lemma provides a bound on (5* by combining Lemma |3] with an analysis of the 
distribution of |W(ti)|. The crucial intuition here is that sets W{el) are strongly 
correlated and hence |VF'(u)| tends to be much smaller than X]j=i 

Lemma 5. SI < 5' := ^MH, + (1 - i) Eg>i j^MH,+,. 

Proof. Let us assume that v has degree d{v) > 3 and that i < d{v), the other 
cases being analogous and simpler. Recall that we need to delete all the edges in 
W^{v) in order to delete ej, . . . , e^, and this takes time M H^^if^y^^ in expectation 
by LemmalU Let us study the distribution of |W^'(w)|. Consider the subtree T' of 
^bin given by (the edges corresponding to) el, ... ^e^ plus their sibling (possibly 
dummy) edges. Observe that, by our assumption on the structure of 5'^j„, this 
tree has i + 1 leaves and height i. We expand T' by appending to each leaf v of 
T' the only path of unmarked edges from v down to some leaf £{v) of S^^^: let 
C" be the resulting tree (with i -I- 1 leaves) . Each edge of T* with both endpoints 
in (the leaves of) C" is a witness edge in W^{v). The number of these edges is 
at most i since the witness tree is acyclic: assume pessimistically that they are 
exactly i. Let r' be the root of T', and s' be the (only) leaf of T' such that the 
edges on the path from r' to s' are unmarked. Let also d' be the only leaf of T' 
which is not the endpoint of any el, j < i {d' is defined since i < d{v)). Observe 
that Pr[s' = d'] ~ 1/2* since this event happens only if the i edges along the 
path from r' to d' are unmarked. When s' ^ d' , there are st most \E{P') \ + 1 
extra edges in W'^{v), where P' is a maximal path of unmarked edges starting 
from r' and going to the root of S^^^. If hi^^ is the maximum value of \E{P')\, 

^° This is irrelevant for [3j , but it is useful in the proof of Lemma [5] 
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then Pr[\E{P')\ = q] = i/2™'i{9+i, for q £ [0,/i,,„]. Altogether: 

9>1 9=0 

The reader may check that the above satisfy Properties (a) and (b) since 
M (i + 1)2^+1 2'+i^^2'?(g + i + l) 2'+!*-^ 2" 

^ ' q>l ' q>l 

is a positive decreasing function of i. Theorem [T] immediately follows from ^ 
and LemmaElsince ^ ^ + 3iJ2q>i " ^ ^ "x) ^ ln4 - |. 
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Figure 1 Decomposition of a tree into components of maximum degree A = 3 
as in Lemma [T] (Left) Part of the edges are labelled with their weight. Squares 
denote terminals. The chosen root is /. The only split node is v. The correspond- 
ing sets of children Vi and V2, for A' — 2, are indicated by the gray ovals. Bold 
edges denote the path Pm{i) — Pi associated with Vi. (Middle) The resulting 
two components. The gray dashed arrow illustrates the charging for the new copy 
of node u, and the black dashed arrow the charging of the nodes in Pi — {v} 
(dashed oval on the right) to the nodes in P2 — {v} (dashed oval on the left). 
(Right) The corresponding component graph. 




Appendix 

Proof. (Lemma\^ Consider the optimal fractional solution x* . We define a fea- 
sible fractional solution x' where x'q ^ = Q for \Q\ > k. Initially x' = 0. For any 
Xq^^, apply the Decomposition Theorem [2] to Q, hence obtaining a collection of 
/c-components Ci, . . . ,Ch- Direct the edges in the component graph towards s, 
so as to identify a sink node Si for each d. For each i, increase aj^^y^j-,.^ ^, by 
Xq j . For a constant k large enough, x' costs at most 1 + e times more than x* . 
Consequently, in order to compute a 1 -I- £ approximate solution, it is sufficient 
to consider the pairs {Q,s) with \Q\ < k, which are polynomially many. The 
number of constraints remains exponential, however the separation problem can 
be solved in polynomial time by the same reduction to MinCut as in [3] . □ 

Proof. (Lemma\Bi) Each iteration takes polynomial time. At any given iteration 
if there is no Steiner tree of zero-cost edges, there exists some terminal r' ^ r 
such that '}Z(^Q,s):seQ(iR.s^r' .r'^Q^%,s ^ 1 ^-^d Pq > for all the considered Q. 
Since w.l.o.g. x*q ^ <1 and hence 1 < X)(q s)^*q s — with probability at 

least l/n'^'^^'> in the current iteration we set to zero the cost of some edge. The 
claim follows. □ 
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Figure 2 Decomposition of a component as in Theorem[21 (Left) A component 
C after contracting nodes of degree 2 other than the root r. Squares denote 
terminals and black nodes are marked in the case g = 1. Dashed lines suggest 
the partition of C into edge-disjoint subtrees. Bold edges indicate the path P{v). 
(Middle) The resulting set of components Ct- regular edges indicate the subtree 
T associated to Ct, and bold edges the paths P{w) associated to the leaves w of 
T. There are two components containing u: the left one because the left child of 
u is along the path P{v) of marked node v, and the right one because u belongs 
to the subtree T of v. (Right) The corresponding component graph. 




Figure 3 (Left) A Stciner tree S*. Squares denote terminals. Edges are labelled 
with their costs. (Middle) The corresponding binary tree S^j^. Bold edges are 
marked. Gray edges define the witness tree T* . The witness sets for the edges of 
cost 8 and 5 are {cd} and {bc,cd,cf}, respectively. Note that these sets have a 
non-empty intersection. (Right) Black edges denote the subtree C' associated 
with the two most expensive edges incident to u, of weight 8 and 5: regular 
edges denote T' and bold edges the paths of unmarked edges from the leaves 
of T' to terminals. The picture also shows the nodes r' , s' , and d' of T'. The 
corresponding witness set is W^{v) = {6c, cd, c/}. Edges bc,cd £ W^{v) have 
both endpoints among the leaves of C. In the example d' ^ s' and P' has length 
(since the edge from r' to its parent is marked): this corresponds to one extra 
edge cf G W^{v). 
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